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A RIGID BODY DYNAMICS DERIVED EROM A CLASS OE 
EXTENDED GAUDIN MODELS : AN INTEGRABLE DISCRETIZATION 

FABIO MUSSO, MATTEO PETRERA, ORLANDO RAGNISCO, AND GIOVANNI SATTA 


Abstract. We consider a hierarchy of classical Liouville completely integrable models sharing the same 
(linear) r-matrix structure obtained through an N-th jet—extension of su(2) rational Gaudin models. The 
main goal of the present paper is the study of the integrable model corresponding to A" = 3, since the case 
N = 2 has been considered by the authors in separate papers, both in the one—body case (Lagrange top) 
and in the n-body one (Lagrange chain). We now obtain a rigid body associated with a Lie-Poisson algebra 
which is an extension of the Lie-Poisson structure for the two-field top, thus breaking its semidirect product 
structure. In the second part of the paper we construct an integrable discretization of a suitable continuous 
Hamiltonian flow for the system. The map is constructed following the theory of Backlund transformations 
for finite-dimensional integrable systems developed by V.B. Kuznetsov and E.K. Sklyanin. 


1. Introduction 


In ^ we have considered Liouville completely integrable Hamiltonian systems with N degrees of freedom 
obtained through the A-th jet-extension of Gaudin models. 

Such procedure allows one to construct a hierarchy of integrable models sharing the same (linear) r-matrix 
structure, whose first elements are, in the su(2) case P E] : 

• A = 1: an Euler top associated with the Lie-Poisson algebra su(2) PIllEl. 

• A = 2: a Lagrange top associated with the Lie-Poisson algebra e(3) = su(2) 0s ppE]. 


Furthermore, a direct sum procedure allows to build long-range chains of n interacting bodies, with rational, 
trigonometric and elliptic i—matrices PP. For instance: 

• A = 1: a Gaudin model associated with 0tisu,(2) PP. 

• A = 2: a Lagrange chain associated with 0"^^ ei(3) PP. 

Moreover this construction can be generalized to any finite-dimensional simple Lie algebra instead of su(2) 

in 


The aim of the present paper is to investigate the one-body system corresponding to A = 3, whose 
underlying algebra still includes 5u(2) and as proper subalgebras but is no longer a semidirect sum of 
subalgebras. This fact suggests one to see the A = 3 case as a slightly generalized version of the Lagrange 
top. The corresponding n-boby system will be considered in a separate paper. 

We have here to mention that, up to our knowledge, the system considered in this paper has been 
introduced, in a different framework, by J.L. Thiffeault and P.J. Morrison in P and it is called the twisted 
Lagrange top. They study this model in the spirit of the dynamical systems theory, so that they do not 
use a Lax pair and an r-matrix approach, as we do in the present work. In P they obtain this new kind 
of integrable top adding a cocycle to the Lie-Poisson structure for the two-field top P, thus breaking its 
semidirect product structure. We remark that in P the so-called twisted top remains a mathematical 
construction without a physical interpretation. 

Later on, O. Vivolo, in nni, constructs a Lax matrix for such system, called here generalized Lagrange 
top. The integrability is proven by direct inspection since an i—matrix approach is not used, and the 
author focuses his attention on the study of the spectral curve of the system through the algebraic-geometry 
machinery. The main goal in nni is the proof that the generalized Lagrange top has monodromy, as well as 
the standard Lagrange top, so that it does not admit global action-angle variables. 

The main purpose of the present paper is the study of these generalized Lagrange tops using the r- 
matrix structure inherited from Gaudin models p. We obtain complete integrability for a large hierarchy 
of integrable systems by providing a Lax pair and a linear i—matrix algebra. In the second part of the paper 
we obtain an integrable discretization of a suitable continuous Hamiltonian flow for the system. The map 
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is constructed following the theory of Backlund transformations for finite-dimensional integrable systems 
developed by V.B. Kuznetsov and E.K. Sklyanin in the papers [niini 

We would like here to remark that our approach allows one to consider also a generic finite-dimensional 
simple Lie algebra instead of su(2) and a natural generalization to a n-body system with rational, trigono¬ 
metric and elliptic dependences on the spectral parameter 


2. Jet-EXTENSIONS OF su(2) RATIONAL GAUDIN MODELS 
Let us consider one-body rational Lax matrices of the following form: 


( 2 . 0 . 1 ) 


^(A) = ^ ^ 


Q!=l 


N-1 




2=0 


/o(A) = i 


where 6“ G K, a = 1, 2,3 and 


M 0 i 

i 0 


^2 = 


0 1 

-1 0 


1x3 = 


1 / i 0 
0 -i 


The 3N coordinate functions yf, a = 1,2,3, satisfy the following Lie-Poisson algebra: 


( 2 . 0 . 2 ) 




■ {y^^yj} 


^ / e“^7 yl+j i + j<N, 

0 i + j > N. 

Here is the skew-symmetric tensor with = 1. Let us notice that N is exactly the order of the 

jet-extension of the Lie-Poisson algebra su(2) It coincides with the number of degrees of freedom of the 
system. 

The functions fi{X) are chosen as 

/ y \ lud 

/o(A), 1=1,...,N-1, 


{q-}fc=l 

where the c^’s are arbitrary complex constants and 


{qJ = {qeN* : qi+ 2 q 2 -\ - \-iqi = i}, |qz| = 


i = l,...,N -1. 


k=l 


In ^ we have shown that the following proposition holds. 

Proposition 1. The Lax matrix IP..0. /I) satisfies the linear i—matrix algebra 
(2.0.3) {'^(A) 01,1® C.{gL)} = [r(A — y), £(A) (8> 1 -I-1 0 -C(y)], 


1 ^ 


A 

Ol — L 

whit respect to the Lie-Poisson algebra Here 1 is the 2x2 identity matrix. 

The Lax matrix (ITint can be written in an equivalent form as a 2 x 2 matrix with elements in the 
negative part of the loop-algebra A“^]: 

i f u{X) v{X) 
w{X) —u{X) 


(2.0.4) 

where 


£(A) = 


N-1 


N-1 


N-1 


i{X) =b^+^ f,{X)yf, u(A) = E/*(A) (j/* w{X) = b^ +ib‘^ +J2 M^) (vl+'^yi)- 


2 = 0 


2=0 


2 = 0 
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It is easy to see that equation (ITTOli is equivalent to the following Lie-Poisson brackets for the rational 
functions u(A), v{X), w{X): 

{u(A), u{n)} = {u(A), u{n)} = {w{X),w{n)} = 0, 


{u{X),v{fj.)} = 
{u{X),w{fi)} = 


X — fj, 


[u(A) -u(/r)], 
[w{X) - w(/r)], 


X — fj, 

{u(A),u;(/x)} = ^^^^[u{X) 


X — fj, 


Recall that the correspondence 


1^3 




^ 2 V 


£ su(2), 


is an isomorphism between 5u(2) and the Lie algebra (K.^, [•,•]), where the Lie bracket [•, •] is realized with 
the wedge product A. Note that 

{^,V) = -2tr(^77) = 2tr(77^*), 
where (■, •) is the scalar product in K^. 

We now focus our attention on an interesting reduction of the Lax matrix If we assume ci = — 1 

and Cfe = 0 for 2 < fc < — 1 we readily get 


(2.0.5) 






N-1 

E 

i=0 


yr 

X^+l 


Let us now fix the following notation: 

y^ = i = o,--,N-i, b = e 

The complete integrability of the hierarchy of systems obtained from the Lax matrices 12.0.511 correspond¬ 
ing to different values of N is established by the following statements. 

Lemma 1. The characteristic curve pl^l : det(/l(A) — = 0 is an hyperelliptic curve of the following 

form: 

(2.0.6) : 4^2 ^ ^b,b) + H{X) + C{X) = 0, 

where 

(2.0.7) 

( 2 . 0 . 8 ) 


N 




2N 


i-2 


2(b,yi_i) +^(yfc,y,_fe_ 2 ) 




iV-1 


c'(A) = E w E {yk,y^-k-2)■ 


X 

i=N-^l k=i — N—l 

Proof: A straightforward computation. 


□ 


Proposition 2. The curve I2.t).h\) provides a set of 2N Poisson-commuting integrals of motion given by 

k-2 N-l 

iJfe = 2(b,yfc_i}+^(yi,yfc_i_2), Ck= ^ (y*,y7v+fe-i-2), fc=l,...,7V, 

2=0 i—k—1 

{H,,Hk} = {Q, Hk] = {Q, Ck] = 0, i,k=l ,..., N. 

The integrals Hk, k = are first integrals of motion. The integrals Ck, k = are linear 

combinations of the Casimir functions of the Lie-Poisson algebra namely 

{Cfe,yf} = 0 Vi/f GflW,Vfc = l,...,iV. 
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Proof: The quantities Hk and Ck, k = 1,N, are immediately obtained through the following formulae: 
Hk = ResA=o Ck = ResA=o A'^+'=-1C(A), k = 1, TV, 

where H{X) and C(A) are given respectively in (12.0. 711 and 112. 0.811 . The function H{X) provides N Poisson- 
commuting first integrals of motion thanks to the r-matrix structure (12. 0.811 . The fact that C'(A) is a 
generating function for the Casimirs of can be proven by a direct computation: 

3 N-l 3 N-l 

{Ck,yj} = {y?y%+k-i-2,yj} = J2 H [y?{y%+k-i-2,yj} + y%+k-i-2{yT,yj} = 

Oi—li—k—1 oc—li—k—1 

3 A^-1 

= 2^ ^ e“^-,2/7+,2/A+fe-.-2+,> 

Oi—1 i—k—1 


for all k = Now, if i + j > N then {Ck,y^}=0 thanks to 12.0.211 . Let us consider 

i+j <N: 




N-l N-l 

'y ^ Vi+j yN+k—i-2+j d" ^ ^ Vi'+j i'—2+j 


a=l 


_i—k—l 


i' — k—l 


where i' = N + k — i — 2. 


□ 


3. The 3-rd jet extension of su(2) rational Gaudin models 


The Lax matrix 


(3.0.9) 




5“ + ^ 
A 


Cl y? y 2 

A2 ^ A2 



is obtained considering TV = 3 in Here 6“ £ M, a = 1, 2,3 plays the role of an external held, taken 

as uniform and constant (in time), Ci,C 2 are real arbitrary constants and the 9 coordinate functions j/“, 
a = 1, 2, 3, T = 0,1, 2 satisfy the Lie-Poisson brackets (innii with TV = 3 (i.e. g^^^), namely: 


(3-0.10) {j/o, 1 / 2 } = 2/1 > 

{y?, 2/f } = yl {//m 2/f } = 0, { 2 / 2 “, 2/f } = 0. 

Let us notice that 

g<^^) = 5 u( 2) 0s 0, dim(y = 6, 

where Q, although including the abelian proper subalgebra K.^ spanned by y 2 , doesn’t have a semidirect 
structure. In it is shown that our algebra g^^^ can be obtained adding a cocycle to e(3, 2) = su(2) 0s 
(R^ 0 R^). We recall that it is possible to use the Lie-Poisson algebra e(3,2) to describe a two-held top [^. 
The Lie-Poisson bracket between two functions /, g : g*^^^ ^ R reads 


{/(yo,yi,y 2 ), 5 (yo,yi,y 2 )} = (yo, Vy^ /a Vy^ g) 0 (yi,Vyo/ a Vy^^T Vy^ /a Vy^g) 0 

+ (y2, Vyo / A Vy, g 0 Vy, / A Vy^ g) 0 (y 2 , Vyi / A Vyi g), 

where V is a gradient with respect to its subscript. The non semidirect structure lies just in the last term 
of the above equation: this term is absent in the e(3, 2) algebra. 

The Lie-Poisson algebra g*-^^ has 3 Casimir functions: 

= (yo,y2) 0 ^(yi,yi}, = (yi,y2), = (y2,y2), 

which differ from the e(3, 2) Casimirs just for the presence of (yo,y 2 ) in 

Specializing Lemma 1 to the TV = 3 case we immediately have the statement 


Proposition 3. The characteristic curve P^^^ : det(£(A) — /r 1) = 0 provides a set of 6 Poisson-commuting 
integrals of motion: 


r(3) : 4g^ + {h,h) + ^ 

A 


Gi 

^ A2 A3 A4 


A5 



(3.0.11) 
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where 

Hi = 2 (b,yo), 

H 2 = (yo,yo) - 2ci (b,yi+C 2 y 2 ), 

H 3 = 2 ci (ci (b,y 2 ) - (yo,yi+C 2 y 2 )), 
Cl = cf (ci + 2 C 2 + 2 , 

C 2 = -2 c? (c 2 C'( 3 )+C'(^)) , 

C3 = 


4. A SPECIAL REDUCTION OF THE 3-RD JET EXTENSION 


Let US now consider the Lax matrix 12.0.511 with TV = 3 and b = (0, 0, b) G To simplify the notations 
we prefer to rename the 9 coordinate functions yf, a = 1, 2, 3, T = 0,1, 2 as j/q = yf = y^ = z°‘, 
a = 1,2,3. Thus we obtain the following Lax matrix: 

3 

(4.0.12) /:(A) = + 


y 


X 


X 




a=l -■ 

Notice that 14.0.1211 is the extension to a third-order pole (i.e. three degrees of freedom) of the Lagrange 
top Lax matrix HE]. Thus, we may expect a generalization of the Lagrange system. 

First, let us notice that the characteristic curve T^^^ : det(£(A) — y 1) = 0 13.0.1111 provides the following 
integrals of motion: 


(4.0.13) 


iLi= 2 (b,y), 7 L 2 = (y,y)+ 2 (b,x}, iLg = 2 ((b, z) + (y, x)), 


(4.0.14) Ci = 2C«=2(y,z) + (x,x}, C 2 = 2 = 2 (x, z), C 3 = = (z, z). 

As in the Lagrange case the third component of y and the euclidean norm of z are constants of the motion. 
Looking at the brackets mM, and taking into account that y and z span respectively 5u(2) and we 
may interpret them as the total angular momentum of the system and the vector pointing from a hxed point 
(which we shall take as (0,0,0) G R^) to the centre of mass of an axially symmetric rigid body, namely a 
Lagrange top. Let us remark that y does not coincide with the angular momentum of the top due to the 
presence of the vector x. We think of x, whose norm in not constant, as the position of the moving centre of 
mass of the global system composed by the Lagrange top and a material point, whose position is described 
by X — z. Here we are assuming that both bodies have unitary masses. The link between these two systems 
is given by integrals Ci and C 2 ■ If we think of a canonical realization of the Lie-Poisson algebra g(3) xnrm 
in term of three canonical coordinates and their conjugated momenta we can immediately argue that the 
vector X must depend on momenta, since 2 ^. We come back to this point in more detail in 

section 5. 

If we look at the first integral 

H 2 = (y,y) + 2 (b,x), 

we immediately see that it coincides with the physical Hamiltonian of the Lagrange top, where the vector y 
is the angular momentum of the spinning top and the vector x describes the motion of the centre of mass of 
the the top on the surface |x| = c, with c constant. In the case of the Lagrange top the equations of motion 
(in the rest frame) with respect to the Hamiltonian H 2 can be written in the following form: 

r y = bAx, 

\ X = y A X, 

which indicate that x rotates rigidly. In the following subsection we shall derive the equations of motions 
for the system described by emphasizing the main differences with the Lagrange case. 
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4.1. A Lax representation. First, it is useful to introduce the following complex generators: 


± 1 _i_ • 2 

y =y ±iy 


± 1 I • 2 

X = X ±1X ^ 


= z^±iz^. 


In term of {y^, y^), x^), {z^, z^) the Lie-Poisson algebra reads 

{y^,y^} = Tiy^, {y+,y~} = -2iy^, 


{y^,x^} = {x^y^} = Ti'^ 


{y+,a; } = {x+,y } = -2i; 


{y^z^} = {z^y^} = Ti V, {y+,z-} = {z+,y-} = -2i 


{y^, x^} = {y+, x+} = {y ,x } = 0, 
{y3, z^} = {y+, z+} = {y-,z-} = 0, 


{x“,z^} = {z“,z^ 


{x'^,x } = —2iz^, 
} = 0, a, /3 = ±, 3. 


The Lax matrix 14.0.1211 has the following simple form: 
(4.1.1) 
where 


y X 

^0) = s +J + ^ 


z 

V’ 


^-2 


b 0 
0 -b 


^=2 


y y 

,+ _,,3 


<¥ = - 
2 


. y ~y 

are all su(2) matrices. 

We choose as a physical Hamiltonian of the system the first integral H 2 '. 


Z = - 
2 


(4.1.2) 


^ = i^f 2 = ^(y,y) + (b,x) = -Tr {y^ + 2BX). 


We see that the kinetic term is given by the norm of the total angular momentum (up to a factor 1/2), while 
the potential energy is given by the projection of the total centre of mass vector onto the external field. 

In the case of the Hamiltonian (in~^ the equations of motion are given by 


(4.1.3) 


y = b A X, 
x = yAx + bAz, 
z = y Az, 


y=[B,X], 

4A < X = [y,x] + [B,z], 
z = [y,z]. 


We immediately see that the vector x does not rotate rigidly, though z still does. Obviously, since |x| is 
no longer preserved, the integral C 2 = 2 = 2 (x, z) does not imply that the angle between x and z 

is constant. Using equations we obtain that the evolution equation for the vector pointing from 

(0, 0,0) S to the position of the material point is given by: 

^ (x — z) = y A (x — z) + b A z 

We conclude this section giving a Lax representation for equations of motion (ESJ). The proof is straight¬ 
forward. 

Proposition 4. The Lax representation for equations mi is given by 


C{\) = -C{\) = [C{\),M{X)], 


Ma) = A + |. 


5. A CANONICAL REALIZATION OF THE REDUCED SYSTEM 

As we have shown in the previous section, our model is a Hamiltonian system. Our aim is now to find 
a coordinate transformation on the symplectic leaves that makes the system canonical. We will use three 
Euler angles 9 S [0, 27r), <f> S [0,27r) and ip € [0, tt) with their canonical conjugate momenta pe^p^ amd p^. 
Our canonical description is restricted to the following symplectic leaf: 


(5.0.4) 


0 = 


= {(y,x,z) G ^(1) ^ 0^ (j(2) ^ 0^ (^( 3 ) ^ ^ 


Proposition 5. A canonical realization of the Lie-Poisson algebra restricted to the symplectic leaf O 
is given by: 


y = ( sincppe -I- cotd cos 


COS( 


■Pij,, — cos (ppe + cot 6 sin 


sm( 


. n ^ ---—" T r-yj ' . ^ P'tjjf P4> 

smO smO ^ 

X = y/2p^ (sin if sin (/ — cos 9 cos if cos (f, — sin if cos (f — cos 9 cos if sin (f, — sin 9 cos if ), 

z = (sin 9 cos (f, sin 9 sin if, cos 9). 
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Proof: It is well-known m that the canonical realization of an angular momentum y and of a constant 
vector z such that |z| = 1 in the euclidean space is given by 

sin(/) 


y = sin (j)pg + cot 9 cos 


coscj) 


, — cos (j)pe + cot 9 sin 


sin0 


PijjtPtf) 


sin 9 

z = (sin 9 cos cj), sin 9 sin (j), cos 9), 

where q = {9,(j),'ip) are the standard Euler angles and p = ipe,P(j>,P'iij) are their canonical conjugated mo¬ 
menta. Let us recall that if f{q,p) and g{q,p) are two arbitrary functions then 


{f{q,p),9{q,p)} = 


df{q,p)dg{q,p) df{q,p)dg{q,p) 


dqi dpi dpi dqi 

Requiring the Lie-Poisson brackets Id. 0. 1011 restricted to the symplectic leaf O 15.0.411 we easily obtain 
X = y/2p^ (sin ijj sin (j) — cos 9 cos t/j cos (j), — sin ijj cos (f) — cos 9 cos sin </>, — sin 9 cos tp). 

□ 

As we mentioned in section 3 we have finally obtained the vector x described in term of canonical 
coordinates 9,(j),ip and the conjugated momentum p^. In particular, we have |xp = > 0. As a 

consequence the Hamiltonian Ti Km describes a non-holonomic dynamics. 

For the sake of completeness we write the first integrals of motion using the above canonical description. 
The physical Hamiltonian Km takes the following form: 


n = ^ 


P%+pI- ‘2P4>P<t> cosd 


2 sin^ 9 


— b \/2p^ sin 9 cos ip- 


We see that the variable (p is cyclic as in the Lagrange case, while explicitly enters in the potential term. 
The remaining first integrals of motion are: 

= 2^2 = P<j>, 


I 2 = -H 3 = y^2p^ [pg sin Ip + {p^ — cos 9) cot 9 cos Ip — p^ sin 9 cos 'ip]+b cos 9. 


6. An integrable discretization through Backlund transformations 


The theory of integrable maps got a boost when Veselov developed a theory of Lagrange correspondences 
[T11[TH| . These maps are symplectic multi-valued transformations which have enough integrals of motion, 
this definition being a proper analog of the classical Liouville integrability. In the main examples, studied 
by him and later by other authors, the integrable maps are constructed as time-discretizations of classical 
integrable models, see, for instance, [I1II3I11. Moreover these correspondences associate to a given solution 
of an integrable system a new solution, a property reminescent of Backlund transformations (BTs) for soliton 
equations. 

In this paper we apply the theory of BT for finite-dimensional integrable systems, developed by V.B. 
Kuznetsov, E.K. Sklyanin and P. Vanhaecke in the relevant papers cniiii Following this approach we 
look at BTs as special Poisson maps. It is possible to find an exhaustive list of the features of these BTs in 
muni. Some of these are the following ones: 

(1) a BT is a Poisson map that preserves the same set of integrals of motion as does the continuous flow 
which it discretizes; 

(2) a BT is given by explicit formulae rather then implicit equations; 

(3) although a BT is multi-valued, it leads to a single-valued map on any level manifold of the integrals 
of motion; 

(4) when one searches for the simplest BT of an integrable system, then one finds a one-dimensional 
family {Bri\p S C} of them. The Backlund parameter g is canonically conjugate to p, i.e. p, = 
—dF/dg with generating function of {Bj^\g G C}. Here p is bound to g by the equation of an 
algebraic curve (dependent on the integrals), which is exactly the characteristic curve that appears 
in the linearization of the integrable system. This property is called spectrality of the BT; 

(5) the explicit nature of a BT makes it purely iterative, so that it is very well suited as symplectic 
integrator for the underlying model. Here the parameter g is an adjustable discrete time step. 
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6.1. One—point BT. First, let us consider the Lax matrix of our reduced system in the form mm . where 
the entries of C{X) are given by 


'a(A) — b —— + + 


X'-' z'-' V X z .^.y^x^z^ 

x'j 2 +^’ ^(A) = —+ ^ + ^, = T + ^ + 

A one-point BT can be defined as the following similarity transform on the Lax matrix /1(A): 

Br, : £{X) ^ M{X;y)C{X)M-^{X,r]) V A G C, 77 G C, 

with some generically non-degenerate 2x2 matrix Ai{X,r]), simply because a BT should preserve the 
spectrum of /1(A). The parameter rj is called Backlund parameter. 

We use “-notations for the updated variables, so that 


Br, : /1(A) ^ £(A) = 


u(X) v(X) 
w{X) —u{X) 


i(A) — ^ + ^ + TT + 


;(A) = 




+ 


X 




X z _ y X z 

X ^ ^ T A^ A 

We shall consider the similarity transformation between /1(A) and /1(A), namely 
(6.1.1) 7 W(A; 7 ;)/ 1 (A) = £(A)M(A; 7 ?) V A G C, 77 g C, 

with the following intertwining matrix Einni: 

A - 77 -l-pg p 

q 1 

Note that the number of zeros of detAd is the number of essential Backlund parameters. Moreover the 
variables p and q are indeterminate dynamical variables. 

Comparing the asymptotics in A —> 00 in both sides of ll^TTll we readily get 

(6.1.3) 


( 6 . 1 . 2 ) 


M{X-y) = 


det A1(A; ij) = X — rj. 


y y 

p = q = 


f = y^. 


2 b' 2b' 

If we want an explicit map from /1(A) to /1(A) we must express q in term of the old variables. To solve this 
problem one can use the spectrality of the BTs [niini Equation ffTT-lTi defines a map Bp parametrized 
by the point P = ( 77 , 77 ) G T^^^. Notice that there are two points on T^^l, P = ( 77 , 77 ) and Q = ( 77 ,- 77 ), 
corresponding to the same 77 and sitting one above the other because of the hyperelliptic involution: 

( 77 , 77) G : det(/l( 77 ) — 77 1) = 0. 

This spectrality property give us the formula Einni 

77 ( 77 ) — 77 w{q) 


(6.1.4) 


q = 


Kv) 


7 ( 77 ) -h 77 ’ 


where 77 and 77 are bounded by the algebraic curve 


1 


-y = 1 \ -h 


Hi 


El 

n3 


Cl , c, 


H-r + 


C 3 


El 

77 if' rj'^ 77 ’ 77 "* 77 '^ 

and the integrals Hi and Ci, i = 1,2,3, are given respectively in (14.0.1311 and (14.0.1411 . Now the equation 
(16.1.Ill gives an integrable Poisson map from /1(A) to /1(A): 

(A - 77 -I- 2pq)[u{X) - g 77 (A)] -l-p 77 ;(A) 


77(A) = 


7 (A) = 


A - 77 

(A — 77 -I- 2pqYv{X) — 2p{X — q + 2pq)u{X) — p'^w{X) 


A -77 


w{X) = 


w{X) + 2 ( 777 (A) — q^v{X) 
A — 77 


Collecting the negative powers of A the above formulae can be rewritten as an explicit map 

Bn ■ (y,x,z) 1—> (y,x,z). 
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given by 


3 = y3 


X +(pq-'n)y 

2 qb, 


y+ = 

5)3 = + py^ — qx~ — q (pq — rj) y~ + 2pqy^ 

3 _ r,2 „.+ 


X = {2pq — rf)x —2px^—p^y'^+pq{pq —r])y —2p^qy^ 
i+ = y+ - ^{pq-y)y- +2qy^, 


(6.1.5) 


3 _ 


z^ — qz -\- 2 pqx^—pq^x -\-px'^pq{pq — p)y -\-ppy'^ -\- 2 ppqy^ ^ 

3 _ ^2 ™— _|_ \ ry-, 9.( rr, ^ _ ,yZ\ n _l_/»0/l»+ _l_0/»O/^/l»3 


= 2qx^ — q^ X + x'^ + rj ^{pq — rj) y + rjy^ + 2riqy 


= {2pq — rf)z — 2p — 2p'^ qx^ + p^ q^ X — p^ x'^ — rj p q{p q — r])y — 

—p^ rjy'^ — 2 rjp'^ qy^. 

The following statement shows how the one-point BT can be written in a symplectic form through a 
generating function. We restrict our BT Brj to a symplectic leaf of the Lie-Poisson structure by fixing values 
of the Casimir functions 

O={(y,x,z)G 0 ( 3 )|^(i)^^^^ C( 2 )= 72 , C'( 3 ) = i}. 

Let us fix the following notation: 

^ = {y\x^zY, ^ = {f,x^zY, 

X~ = {y~ ,x~ ,z~f', x+= (y+,i+,5+)'^. 

Proposition 6. The one-point BT Br^\o can be arranged in the form: 

3 

( 6 . 1 . 6 ) = 


1=1 

3 


(6.1.7) 

1 = 1,2,3, where 

Fr,{x~\x") = 


^* = E{7^^*}V^-+F,(x-|x+), 

1=1 


-I-+, _ y , ,Jy , y^\ (1 + 1772 )^ ,^(12 


( 6 . 1 . 8 ) 

with 


2 b 

1 

'2fc 




Akr]'^ 




k + 1 
k-1 


_, --i__ X X n _, 

' rn _L -y nr^' _ 'n nr'' _L _ f _ _L _ nr 


Z~^ X Z X —px X + 


z” V 2: 


- X 2:"^ - 72 + 


x'^ f x^ n ^ , 

- ' - + ^X+ Z- -P2 


z+ \ z+ 2 


k^ = l + yz z~^. 

Proof: The Casimir functions do not change under the map: 

) ^ ((7(^7 (7(2), (7(3)) = ((7(^7 , (7^3)). 
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The above invariance allows one to exclude 6 variables, expressing y~^,x^, z~^ and y ,x ,z in term of the 
components of the vectors and 


y = — 


1 


z' = 


1 - (z3) 


2 j2 -2 x^ - (zy) 

z 

312 


y = TT 


2 71 - 1 -2^353 _y(i_ (^3)2) 


X = 


z = 


1-(53) 


2 72 - 25 ^ 53 _ - (1 _ (- 3 ) 2 ) 

312 


With the help of Iti.l.dH we can rewrite equations ( 16 . 1 . 51 ) of the map in the following form: 


'hi = 2/" = 


1 n 


y y^ 


2 b 2 /c 3 1 2 


[y{ 


Z X X z 


)-z-z+]^ + [(z-fz+ y+ + Cz+fz- y-] + 


+ ^1 +-1772^(2 x'^ + X 5 +) — (72 + 2 ?77i) 2; 5 +- 77 (x x^ + z y'^ + z'^ y ) + 
+ ^(7|-47i)|, 


^'2 = = 


— [i/+ X + 77 (x + x+ z ) - 5 + 2 + 72 ] , 


^3 = ^" 


y+ z~ 


2 w 


k, 


'hi = r 

’2 


= y 


^2 = 


^3 = 5 " = 


2 k 
z^y~ 

2 w 


[x+ y + 77 (x 5“*" + x“*" z ) — 5“*" z + 72 ] , 


k, 


where = 1 + 77Z z+. It is now easy to check that the function F„(x Ix'*’) Hb.l.^H satisfies equations 

11331), JQ-TI) . 

□ 


Remark 1. It is possible to use the theory of canonical transformations to show that Bri has the spectrality 
property. The spectrality property of a BT means that the two coordinates 77 and /i parametrizing the map 
are conjugated variables, namely 


y = 


dF^ 

' dy ’ 


where F^j is the generating function of the BT. 

In our case, using equations 11331, (i?n~m and 11331, we obtain 


y~ 

77 = 77 ( 77 ) - -^v{y) = - 


9Fr,{x Ix"^) 

dy ’ 


so that the spectrality property holds. 


Remark 2. We have here to remark that the one-parameter map 116.1.51) is a complex transformation, 
namely it is not a physical BT for our system. Following |2] we shall construct, in the next subsection, a 
physical map using two Backlund parameters. 























A RIGID BODY DYNAMICS DERIVED FROM A CLASS OF EXTENDED GAUDIN MODELS 


11 


6.2. Two—point BT. According to HIIHl, we construct a composite map which is a product of the map 

—Ai2)* 

= Bq, o Bp, : T(A) ^ C{X) ^ £(A). 


The two maps are inverse to each other when rji = r]2 and fJ-i = fJ.2- This two-point BT is dehned by the 
following discrete-time Lax equation: 

(6.2.1) M{X-,r]i,r]2) C{X) = C{X) M{X-,r]i,r]2) VA G C, 771,172 GC, 

where the matrix At (A; 771,772) is piTT)] 


( 6 . 2 . 2 ) 


M{X;r]i,r]2) = 


X — rji + st t 

-I-(771 - 772) s X-r] 2 -st 


det At (A; 771,772) = (A - 77i)(A - 772). 

The spectrality property with respect to two fixed points (771,771) G T^^^ and (772,712) G T^^) give 

Hm) - fJ'2 


(6.2.3) 


(6.2.4) t = 


v{m) v{r]2) 

ivi - V 2 ){u{m) + - M2) 


(771 - r]2){u{m) - Mi)(u(M2) + M2) 


(77(771) -I- 771)777(772) - (77(772) - 772)777(771) (77(772) -I- 772)777(771) - (77(771) - 771)777(772)' 

Now we have two Backlund parameters 771,772 G C. It is possible to obtain several equivalent formulae Emu 
for the variables s and t since the points (771,771) and (772,772) belong to the spectral curve i.e. are bound 
by the following relations 


-(2mjT =P{Vj) + v{Vj)w{r]j) = 772(77^) -h 77(77^)77;(?7j), j = 1,2. 

Together with (16. 2. 311 and (EHIl, the formula gives an explicit two-point Poisson integrable map 

from /1(A) to T(A). The map is parametrized by the two points Pi and Q2- 

Obviously, when 771 = 772 (and 771 =712) the map turns into an identity map. As we have shown in [2] the 
two-point BT can be reduced to a real Poisson integrable map if the following condition holds: 


m = 'n 2 = v = ^iv) + i 3(77) e C. 


Therefore, the two-point map leads to a physical BT with two real parameters. 
Let us now introduce the following notation 

^ = (y, X, z)^ G K®, X = (y, X, z)^ G K®. 

A direct computation based on the similarity transform Km shows that 


Proposition 7. The two-point BT Bp,^Q^\rf,=f )2 • ^ ^ ^ is given by 

(6.2.5) A = $(s,t;77)A + Ao(s,t;77), 

with 

/ l3x3 ©3x3 ©3x3 \ 

$(s,t;77)= A{s,t]r]) I 3 X 3 © 3 x 3 , 

\ B{s,t-r]) A{s,t;v) l3x3 / 

where A{sA\ii) B(s,tiv) ®7’e two 3x3 dynamical matrices depending on the Backlund parameter rj and 
the parameters s,t and Xq(sA',v) is a dynamical vector. The matrices I 3 X 3 and © 3 x 3 are respectively the 
3x3 identity matrix and the 3x3 zero matrix. 


The explicit expressions of Xo(s, t; 77), A(s, t; 77), i?(s, t; 77) are rather complicated and they are given in 
Appendix 1. 

Remark 3. Notice that, despite its matrix formulation, the map Km is a non linear transformation. 
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7. Concluding remarks 

We have considered a hierarchy of classical Liouville completely integrable models sharing the same (linear) 
r-matrix structure obtained through an A^-th jet-extension of su(2) rational Gaudin models. The general 
procedure of such extension is presented in 

We have fixed TV = 3 obtaining a rigid body associated to a Lie-Poisson algebra which is an extension 
of the Lie-Poisson structure for the two-held top. We have here to recall that this classical system has 
been introduced in where it is called the twisted Lagrange top, and furtherly investigated in m in the 
algebraic-geometry setting. 

The novelty of our approach is the introduction of an r-matrix formulation for this system. Its knowledge 
enables us to easily hnd a Lax formulation for the equation of motion, as well as a canonical realization in 
terms of Euler angles. Finally, through the approach developed by V.B. Kuznetsov, E.K. Sklyanin and P. 
Vanhaecke mma, we hnd explicit BTs for the system. In Appendix 2 we present a numerical simulation 
of the real reduction of the two-point BT. 

Let us remark that another feature of our approach is the natural possibility of constructing n-body 
integrable chains starting from each A^-th jet-extension, both considering any simple Lie algebra and rational, 
trigonometric and elliptic dependances on the spectral parameter. This was done for the N = 2 case in [?]. 
The case N = 3 will be considered in a separate paper. 


Appendix 1: Explicit expressions of Xo{s,t,r]), A{s,t-,r]), B{s,t;r]) 


As we have shown in section 6.2 the two-point BT can be written in the following form: 

Br, ■■ X I—> X = $(s, t;r])X + Xo{s, t] rf), 

where 

/ 13x3 ©3x3 ©3x3 \ 

13X3 ©3x3 , 

\ B{s,t;il) 13X3 / 

and Xo(s, t; rj) is a dynamical vector depending on ij, s, t just as like the 3x3 matrices A(s, t; ij) and i?(s, t; rj). 

Here we recall that the dynamical variables s and t can be obtained using the spectrality property of the 
BT and are given by formulas 16.2.311 and 16. 2. 411 . namely: 

„ _ ivi - V 2 ){u{r]i) + - ^^ 2 ) 

vijli) ’ (m(i7i)+Mi)w(?? 2 ) - (u(?72)-M2)w(i7i)' 

The explicit expressions of Aro(s, t; rf) is given by 


[Xo{s,t;'n)]i =b{sai- t), 
[Xois,t;r ])]2 = -b{s ai + t). 


[Xo{s,t]v)h = 


s ai 


[Xo{s,t;v)]4 = b +^{r])) +t{a2 -2IH(r?)) 

[Xo{s,t;v )]5 = ife “^(“2 + ^iv)) -t{a2 - 291(77)) 


[Xo(s,t;77)]6 = 2btsai, 

[Xo{s,t;v )]7 = [(3(17) +^{v)f{t + ai) + 8 stmi'n){t-ai)] , 

[Xois,t;v )]8 = -^ - ai) + 8 st93(77)(t + oi)] , 

[Xo{s,t\r ])]9 = 2 [Xo(s,t; 77)]6 93(77), 


where we have introduced the quantities 


Q!i = 23{r]) — st, 
02 = 0:1 — s t. 


and 3(r]),d\{ri) denote respectively the imaginary part and the real part of the Backlund parameter rj. 
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The entries of the dynamical matrix v4(s, t; rf) are given by: 

[ A { s , t \ ri)]ii = {), i = l,2,3, 

[A{s,f,'n)\i2 = ia2 = -[A{s,t;v)]2i, 

[MsA;v)]i 3 = sai-t = -[A{s, f,v)] 3 i, 

[^(s, t; r])]23 = -i (s ai + t) = -[A(s, t; r])]32, 

so that A(s, t; 77) is a skew-symmetric matrix. The matrix B{s, t; 77) has a more complicated form: its entries 
are 

77)]ii = i(l - s^){al - t ^), 

[B{sA;v)]i 2 = - 2a29^(i?)], 

[B(s, t; 77)]i3 = is ai (02 + ^^(77)), 

[B{s, t- 77)]2i = - s2 ai) + 2 02 ^(77)], 

[B{s,f,v)]22 = i(l + s2)(a? -f) 

[B(s, t; r ))]23 = -^sai (02 + ^{v)), 

[B(s, t; 77)]3 i = i [t (q;2 + 2 9^(77)) + s Qfi (02 - 2 91(77))], 

[B{s, t; tj )]32 = i [t (02 + 2 9 ^( 77 )) - s Qfi (02 - 2 91 ( 77 ))], 

[B(s,t; 77)]33 =Astai. 

Appendix 2: Numerics 

In this appendix we present a 3D plot corresponding to the real reduction of the two-point BT (16. 2. 511 . It 
is obtained using a MAPLE 8 program that is a slightly different version of the MATLAB program developed 
by V.B. Kuznetsov in |5]. 

The input parameters are: 

• the intensity of the external field, i.e. b; 

• the Backlund parameter 77 = 93(77) + Here 3(77) is the time-step of the discretization; 

• the number of iteration of the map, N; 

• the initial values of the coordinate functions y^, x^, z^,z^, z^. 

The output is a 3D plot of N + N consequent points (x^ — z^,x^ — z^,x^ — z^) and z^). We remark 

that the vector (x^ — z^,x'^ — x^ — z^) describes the position of the material point, as explained in section 
4, and the vector (z^, z^) is the position of the centre of mass of the spinning top. As expected, the points 

(z^,z^,z^) lie on the sphere = (z,z) = constant, of some radius defined by the initial data. 
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Figure 1. 


(6;lH(7?),a(77);V;y\y2,y^,x\x2,a;^,2\2^,2^) = (l;5,0.1;1000;-2.4,-0.6,-1.2,-2.19,0.89,l-34,l,0,0) 


-3-2-10 1 2 3 



Figure 2. Projection of Figure 1 on the x — y plane 
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